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Entanglement lies at the heart of quantum mechanics—as a fundamental tool for testing its 
deep rift with classical physics|3lHl while also providing a key resource for quantum technolo¬ 
gies such as quantum computatioiJ^ and cryptography]^ In 1987 Greenberger, Horne, and 
Zeilingei® realized that the entanglement of more than two particles implies a non-statistical 
conflict between local realism and quantum mechanics. The resulting predictions were exper¬ 
imentally confirmed by entangling three photons in their polarization.!^ Experimental efforts 
since have singularly focused on increasing the number of particles entangled, while remaining 
in a two-dimensional space for each particle.®® Here we show the experimental generation of 
the first multi-photon entangled state where both—the number of particles and the number 
of dimensions—are greater than two. Interestingly, our state exhibits an asymmetric entan¬ 
glement structure that is only possible when one considers multi-particle entangled states 
in high dimensions.!^ Two photons in our state reside in a three-dimensional space, while 
the third lives in two dimensions. Our method relies on combining two pairs of photons, 
high-dimensionally entangled in their orbital angular momentum, in such a way that informa¬ 
tion about their origin is erased.!^ Additionally, we show how this state enables a new type 
of “layered” quantum cryptographic protocol where two parties share an additional layer of 
secure information over that already shared by all three parties.!^ In addition to their applica¬ 
tion in novel quantum communication protocols, such asymmetric entangled states serve as a 
manifestation of the complex dance of correlations that can exist within quantum mechanics. 

In the recent past, big mysteries in quantum mechanics have been illuminated by taking small steps in 
the right direction. The phenomenon of quantum interference, for example, appeared when one considered 
a single quantum particle with at least two discrete levels, or dimensions. Moving to two particles gave 
us quantum entanglementfH and Bell’s inequalities,® which allowed the conflict between quantum mechanics 
and local realism to be tested in a statistical manner.!!^ Increasing the number of entangled particles to 
three, while seemingly a simple step, provided the first “all-or-nothing” test of local realism.® Alongside 
this, increasing the dimensions of a single quantum particle from two to three provided the first clear test of 
quantum contextuality.!^ History would dictate that increasing both the number of quantum particles and 
the number of dimensions in concert will lead to further interesting and fundamental quantum phenomena. 
In this letter, we discuss one such phenomenon that arises when we consider three-particle entangled states 
of dimension greater than two—namely, asymmetric multi-particle entanglement. 

Two-dimensional entangled states have been studied in great detail due to their natural application 
in quantum information and computation. However, the amount of information carried by a photon is 
potentially enormous, and the ability to harness this information leads to quantum communication systems 
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with record capacities and unprecedented levels of security!!^ A natural space for exploring large dimensions 
in a quantum system is a photon’s spatial degree of freedom.^ The orbital angular momentum (0AM) 
of a photon is a spatial property that provides a discrete and unbounded state spaceJ^ The phase of an 
OAM-carrying photon winds azimuthally from 0 to 27ri around the axis. The number of “twists” in the 
wavefront are given by the quantum number £ and dictate the photon state dimension. Recent experiments 
have shown the entanglement of two photons in up to 100 x 100 dimensions in their spatial modesEHEHl The 
dimensionality of two-photon entangled states is given by a single number, the Schmidt number,!^ which 
is the rank of the reduced single particle density matrix. This number represents the minimum number of 
levels one needs to faithfully represent the state and its correlations in any local basis. When one considers 
three entangled photons of dimensionality greater than two, the question of how many levels per photon are 
involved has three answer s.^^^ In our experiment, we aim to create a state which has the form 

1^)332 ~ [ I^)a l^)c + I^)a l^)c + I^)a I^)b l^)c ] • (^) 

Notice that the first two photons, A and 5, live in a three-dimensional space, while the third photon, C, lives 
in a two-dimensional space. The dimensionality of this state is given by a vector of three numbers (3, 3, 2). 
Specifically, these are the ordered ranks of the single particle reductions of the state density operator: 

rank(p^) = 3, rank(p 5 ) = 3, rank(pc') = 2, (2) 

where pi = Tr^ 1^)332 (T| 332 is the state of system i G (A^B^C). Only certain combinations of these 
three ranks are allowed, leading to a rich structure of possible asymmetric high-dimensional multipartite 
entangled states.l^ In addition to being of fundamental interest, such states enable a novel “layered” quantum 
communication protocol that allows the sharing of secure information asymmetrically amongst multiple 
parties.l^ We discuss this protocol in the context of our state later in the text. 

The entanglement of three photons was first achieved by combining two pairs of polarization-entangled 
photons in such a way that it became impossible, even in principle, to know which pair one of the detected 
photons belonged toUillHHl The workhorse of such experiments is the polarizing beam-splitter (PBS), which 
was designed to separate a light beam into its horizontal and vertical polarization components. Interestingly, 
a PBS can also be used to mix the polarizations of two input photons from independent polarization- 
entangled pairs in such a manner that an output photon contains polarization components from both input 
photons. In this manner, information about their origin is erased, producing a four-photon, two-dimensional 
GHZ state.l^ The largest such state created thus far used seven BBSs to entangle eight photons in their 
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Figure 1: Schematic of the experiment. A blue laser pulse creates two pairs of OAM-entangled photons in two nonlinear 
crystals (NLl and NL2). Both pairs are spatially separated via polarizing beam splitters (PBS). A moving trombone prism 
ensures that a photon from each pair arrives simultaneously at the OAM beam splitter. This device contains two dove-prisms 
(DPI and DP2). DPI reflects an incoming photon that inverts its spiral phase front, while DP2 both inverts and rotates it 
by 180°. This results in interference that depends on the parity (odd or even) of the spatial mode. Considered as a whole, 
the OAM beam splitter reflects odd components of OAM and transmits even ones, mixing the OAM components of two input 
photons. A coincidence at detectors B and C can only arise when both detected photons have the same mode parity. This 
projects photons at detectors A, B, and C into an asymmetric three-photon entangled state that is entangled in 3 x 3 x 2 
dimensions of its OAM. The detection of a trigger photon at detector T heralds the presence of this state (see extended data 
flgure 1 for an in-depth schematic). 


polarizationP In order to manipulate the high-dimensional space of OAM, one would need a device akin to 
the PBS, but operating on a photon’s spatial wavefunction as opposed to its polarization. 

In our experiment we choose to use a simple Mach-Zehnder interferometer that was designed to sort 
photons based on the parity of their OAM quantum numberAs shown in Fig. [T] this interferometer has 
dove prisms (DPI and DP2) in both arms, rotated by 90° with respect to one another. DPI reflects an 
incoming photon that inverts its spiral phase front, while DP2 both inverts and rotates it by 180°. Thus, 
each input photon interferes with a rotated version of itself, leading to constructive or destructive interference 
based on whether its spatial mode has odd or even parity. We use this device as a two-input, two-output 
“OAM beam splitter” that reflects photons with odd OAM values and transmits even ones. In this manner, 
an outgoing photon contains a superposition of odd and even values of OAM from two independent input 
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photons. 


Our procedure for generating the high-dimensional tripartite state uses two independent entangled photon 
pairs created in two non-linear crystals (NLl and NL2). Suppose that the two pairs are in the state 

\'^)abcd — \~'^)b + |0)a |0)b + I"^)a I^)b ) 

1 "^)D + l^)c + l"^)c I^)d ) ’ (^) 

which is a tensor product of two OAM-entangled photons pairs (A, B and C, D) of dimension d = 3 each 

(0AM quantum number £ G { — 1,0,1}). At this stage, mixing these two states results in nine possible 
four-photon amplitudes. Photons B and C are then incident on the two inputs of the 0AM beam splitter. 
Since this device reflects odd values of 0AM and transmits even ones, a coincidence detection between 
the two outputs can only arise when photons B and C are both carrying either odd or even values of 
0AM. This projects the state from Eq. Q onto a subspace spanned by the five terms (|l)y^ \~'^)b I^)c 

(|1)a I"1)b I"^)c I^)d)’ (I^)a I^)b I^)c 1 ^) 1 ))’ (I"^)a I^)b I^)c I"^)d)’ (I"^)a I^)b I"^)c I^)d)- Photon D is 

then measured in a superposition state given by ^ mode-projection carried 

out by a spatial light modulator (SLM) and a single-photon detector (T). This acts as a trigger for the 
three-photon entangled state: 

l^)exp = ^ [|1 )a |-1>b |1)c + |0U |0)b |0)c + |-1>^ |1)b |1)c]- (4) 

Note that this state has the same form as equation except photons A and B are in a three-dimensional 

space given by the 0AM quantum numbers (-1, 0,1), while photon (7 is in a two-dimensional space given by 
(0,1). In order to ensure that these three terms are in a coherent superposition (as opposed to an incoherent 
mixture), we perform superposition measurements in a two-dimensional subspace spanned by the second and 
third terms in equationMeasuring photon A in a superposition state given by = :^(|0 )a “ I“^)a) 

projects the three-photon state into: 

T (|p)^' |M)y + |M)^' |p)y), (5) 

where = ~^iA)B/c + |4)b/c) \^)b/c ~ ^(I*^)b/c “ |4)b/c)- terms |P)g |P)c> and 

iM)y |M)y are missing because of two-photon destructive interference at the 0AM beam splitter, which 
only occurs when photons B and C are indistinguishable.!^ Note that this requires both two-photon and 
one-photon interference (for two different photons) to occur at the 0AM beam splitter, which proved to be 
an experimental challenge. The use of narrowband interference filters (IF2) before detectors B and C blurs 
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out the temporal correlations of the photons with their entangled partners to a certain extent, ensuring that 
they cannot be distinguished based on their time of arrival!^ We look for destructive interference between 
photons B and C by projecting them into states |M)^^ and |M)^^ , and scanning the trombone prism shown 
in Fig. At the zero-delay position, both photons arrive at the first beam splitter (BSl) at precisely the 
same time and interfere destructively, leading to a drop in coincidence counts between detectors B and C, 
conditioned on photons A and D being detected in states and |P)^'^ respectively. The characteristic 

dip in Fig. [^has a width of 473/im that matches what we expect from the bandwidth of our filters and the 
thickness of our crystals (see Methods section for details). 

In order to verify that our state is indeed genuinely multipartite entangled in a 3 x 3 x 2 dimensional 
space, we have to prove that it cannot be decomposed into entangled states of a smaller dimensionality 
struct ure}Hll2II such as a (322)-state or a two-dimensional GHZ/(222)-state. As shown in Fig. [^a), these 
states take on a “russian doll” structure of concentric subsets embedded in the set of (332)-states, and can 
be ruled out as follows. First, we find the best achievable overlap F^ax of a (322)-state a with an ideal 
(332)-state |T). In the Methods section, we show that Fmax = 2/3. Next, we calculate the overlap Fexp of 
our experimentally generated state pexp = |^)exp (^lexp with an ideal (332)-state. If Fexp > Fmax, our state 
is certified to be entangled in 3 x 3 x 2 dimensions. In order to calculate Fexp, it is sufficient to measure the 



Figure 2: Three-photon coherent superposition. Experimental data showing the result of measuring photons A, and C 
in states and |M)^^ respectively, conditioned on measuring photon D in state |P}^ ^ • The drop in coincidence 

counts at the position of zero delay results from interference between photons B and C. This indicates that the 3-photon state 
from equation ^ is in a coherent superposition. A Gaussian fit is calculated based on the spectral shape of our narrowband 
filters and has a full-width at half-maximum of 473)Um and a visibility of 63.5% (see methods section for details). 
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18 diagonal elements of Pexp, as well as its 3 unique off-diagonal elements that are expected to be non-zero. 
The 18 diagonal elements are measured directly through projective measurements of 27 minutes each. The 
3 off-diagonal elements are measured via a series of 144 projective measurements (see Methods section for 
details). These elements are plotted in Fig. |^b), and numerical values are provided in extended data table 
1. From these, we calculate an experimental fidelity Fexp = 0.801 ± 0.018, which is above the bound of 
Fmax = 0.667 by 7 standard deviations. This certifies that our three-photon state is indeed entangled in 
3x3x2 dimensions of its orbital angular momentum. 



Figure 3: Witnessing genuine multipartite entanglement in high dimensions, (a) In order to verify that our three- 
photon state is entangled in 3 X 3 X 2 dimensions, we have to show that it cannot be decomposed into entangled states of a 
smaller dimensionality structure. First, we calculate the best achievable overlap of a (322)-state a with an ideal target (332)- 
state 1^} to be Fmax = 2/3. Next, we calculate the overlap Fexp of our experimentally generated state pexp with the target 
state 1^}. (b) The 18 diagonal and 3 unique off-diagonal elements of pexp that are measured in order to calculate a value of 
Fexp = 0.801 ± 0.018 (elements not measured are filled in with lines). This is above the bound of Fmax = 0.667 by 7 standard 
deviations and verifies that the generated state is genuinely multipartite entanglement in 3 X 3 x 2 dimensions of its orbital 
angular momentum. Fgxp does not reach its maximal value of 1 because the state superstition is not perfectly coherent. 
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Figure 4: A layered quantum communication protocol, (a) When Alice, Bob, and Carol share an asymmetric (3, 3, 2) 
entangled state, all three have access to an entangled bit in the (0,1) basis. Additionally, whenever Carol has a 1, Alice and 
Bob share an entangled bit in the (1,2) basis, (b) After sufficient data is collected, these bits can be used as a one-time pad to 
encrypt a secret message shared among all three parties. In addition, Alice and Bob can share a second layer of information 
unknown to Carol. Security is verified by checking for the presence of (3, 3, 2) entanglement in a randomly selected subset of 
photons. 


Finally, this high-dimensional three-photon entangled state enables a new type of “layered” quantum cryp¬ 
tographic protocol where layers of information can be shared asymmetrically with unconditional security.l^ 
Consider the communication scheme depicted in Fig. If Alice, Bob, and Carol share a (332)-entangled 
state, all three parties have one entangled bit in the (0,1) basis. However, Alice and Bob have access to 
an additional entangled bit in the (1,2) basis whenever Carol has a 1. This state allows all three parties 
to generate a secure string of random bits that can be used as a key. Two-thirds of the time, however, 
Alice and Bob are able to generate a second secure key. Thus, all three parties can securely share a secret 
message, with Alice and Bob sharing an additional layer of information that is unknown to Carol. Security 
can be tested by verifying the presence of (332)-entanglement in a randomly selected subset of photons using 
the entanglement witness from above. The implementation of such a protocol would require multi-outcome 
measurements of 0AM, techniques for which are quite mature today.^^ 

We have created a three-photon entangled state and verified that it is high-dimensionally genuinely 
multipartite entangled in its orbital angular momentum. This is the first such entangled state where both. 
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the number of particles and the number of dimensions, are greater than two. In addition, our generated state 
displays a unique asymmetry in its entanglement structure that has not been experimentally observed before. 
This asymmetric entanglement has enabled us to develop a new layered quantum communication protocol 
where different layers of information can be shared securely amongst multiple parties. Our experimental 
method can be modified to create a vast array of even more complex entangled stateP^^ —of which the state 
above is just one example. Such asymmetric entangled states constitute a new direction in experimental 
studies of entanglement, and will allow the development of complex, multi-level quantum networks in the 
future. 
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Methods 

Down-conversion sources. As shown in extended data Fig. 1, a femtosecond pulsed TiiSapphire laser 
(Coherent Chameleon Ultra II) with a center wavelength of 808nm is focused by lens LI into a ^d-Barium 
Borate crystal (BBC) to generate blue pump pulses at 404nm via the process of second-harmonic generation. 
The BBO crystal is incrementally moved by a motor (M) in order to avoid damage to it. The generated 
pump pulses (average power of 820mW) are re-collimated by lens L2 and separated from the 808nm light by 
a dichroic mirror (DM1). They are then focused with lens L3 onto two periodically poled potassium titanyl 
phosphate (ppKTP) crystals (NLl and NL2) with dimension 1mm x 2mm x 1mm and poling period 9.55/im 
for type-II phase matching. Two independent photon pairs entangled in their orbital angular momentum 
(0AM) are generated in the two crystals via spontaneous parametric down-conversion of the same pump 
pulse. The crystals are spatially oriented so that down-conversion occurs when the pump pulses are vertically 
polarized. In order to conform to the phase-matching conditions for producing wavelength-degenerate pho¬ 
ton pairs at 808nm, NLl (NL2) is heated to 119.1°C (121.9°C). A 4/ imaging system composed of two lenses 
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(L4) between between the two crystals ensures that the waist of the pump pulse is approximately 240/im at 
both crystals. The pump waist is chosen such that generation rate of / = 0 photons is approximately twice 
that of the / = ±1 photons in the two OAM-entangled photon pairs. Note that the 0AM state distributions 
in equation © are flat, while our experimentally generated distributions are not. We account for this by 
unbalancing the trigger photon (D) state superposition to = 0.5110)^ +0.86 |-1)/), which surprisingly 

lowers our average four-fold count rates only by 6%. Two dichroic mirrors (DM2) separate the 404nm pump 
pulses from the entangled photon pairs. Due to the high pump power, the Kerr-lensing effect cannot be 
neglected and adds an equivalent lens with Im focal length into the optical path. This is compensated 
by carefully aligning the 4/ imaging system (L4). To avoid gray-tracking due to the high pump power, 
the ppKTP crystals are also moved incrementally in steps of 200/im with motors (M). Movement of the 
BBO and ppKTP crystals shows slight variation in the count rates, which is compensated for by averaging 
over several periods of movement within each measurement interval. Two polarizing beam splitters (PBS) 
deterministically separate both photon pairs such that all four photons are propagating in different directions. 

OAM beam splitter. This is an interferometric two-input two-output device that mixes input photons 
based on the parity of their OAM mode. In Fig.[^in the text, the OAM beam splitter is depicted for simplic¬ 
ity in its original desigrP^ as a Mach-Zehnder interferometer. However, our experiment requires long-term 
interferometric stability, which is not easily attainable with a Mach-Zehnder interferometer. For this reason, 
we implement the OAM beam splitter in a zero-area, double Sagnac configuration. Highlighted in green in 
extended data Fig. this configuration takes a standard Sagnac loop (which is one-input and one-output) 
and shifts it laterally, creating two counter-propagating Sagnac loops side-by-side that meet back up at the 
beam splitter. This allows us to install a dove prism (DPI and DP2) in each loop. The interferometer is 
further implemented in a zero-area configuration that is rotationally invariant, minimizing the rotational 
dependence of the entire interferometer on the rotation of the earth. Furthermore, a piezo controlled mirror 
(P) is scanned every hour to optimize the alignment of the interferometer and ensure it is sorting OAM 
modes correctly. In addition to controlling one of the mirrors actively, the interferometer is enclosed in a 
plastic box that reduces the air flow through the device and stabilizes the temperature to within ±0.02°C. 
All of these efforts lead to a very long stability (> 80 hours) and at the same time achieve a very high sorting 
efficiency of even and odd OAM modes (99:1). 

Interference of independent photons. The interference of two photons at the OAM beam splitter 
crucially depends on the indistinguishability of the two particles involved. Since the presence of the interfering 
photons is indicated by the trigger detection events of their respective partner photons, we effectively observe 
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four-photon coincidence events. The interference visibility strongly depends on the physical properties of 
these four photons, as well as the process by which they are created. This particular phenomenon has been 
studied in great detaiP^ and the theoretically expected visibility is calculated to have the form 


y' = 


^ a/^T + ^pV^s + 

crpy/crl + + cr|> 


(6) 


Here, ap, and crp are the gaussian spectral widths of the pump, signal, and trigger photons respectively 
(in units of frequency), as and crp are determined by the widths of narrowband interference filters IF2 and 
IFl, respectively, tj has units of time and refers to the relative timing jitter between the two crystals, which 
arises due to the group velocity mismatch between the 404nm pump pulse and the 808nm down-converted 
photons in the ppKTP crystals. Additionally, due to the finite pulse width of the 404nm pump beam, a 
second relative timing jitter between the two crystals is also present. In our experiment, the sum of these 
two timing jitters is approximately Ips.^^ 

In addition to spectral distinguishability, the spatial and temporal distinguishability of the two photons 
both play an important role in the interference visibility. The moving trombone system of mirrors (TB) is 
used to find the position of zero-delay when photons from both crystals arrive at the 0AM beam splitter 
at the same exact time. The spatial-mode distinguishability of the two input photons at the 0AM beam 
splitter has a further detrimental effect on the total interferometric visibility. We use a 4/ imaging system 
of lenses (L6) in order to compensate for the extra propagation that the photon from NLl undergoes with 
respect to the photon from NL2. Nonetheless, the mode overlap of the two photons at the 0AM splitter is 
not perfect, and modifies the two-photon interference visibility in the following way: 


r?V' 

V= - - -. (7) 

i + w(i-7?2) o; 

Here, r] = t^oam x ^sp is the product of the sorting efficiency t^oam of tho 0AM beam splitter and the spatial 
mode overlap r/sp of the two input photons. Inserting our measured/estimated experimental parameters 
{ap = 3.67THz, ap = 588GHz, as = 184GHz, rj = Ips, 770 AM = 0-99 and 77 SP = 0-9) into the above 
formula yields a two-photon interference visibility of V = 0.64, which is very close to the experimental 
visibility observed in Fig. The visibility crucially depends on the spectral width as of the two interfering 
photons. A smaller as greatly improves the interference visibility at the cost of lowering the four-photon 
count rates. It is therefore important to find an optimal value that provides both, sufficiently high visibility 
and practical count rates. In our experiment, we found this value to be as = 184GHz. Assuming a Gaussian 
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spectral distribution, the theoretical full-width at half-maximum (FWHM) of the two-photon interference 
dip is then given by 


-^dip — \ 


21n2 


TT \/ cr 


-2 


+ r 


2 ' 


( 8 ) 


Inserting the corresponding values from above leads to a expected FWHM dip-width of 624/im, which is 
close to our experimentally observed value of 473/im in Fig. 


Single Photon Projective Measurements. Two additional 4/ imaging systems (L6) image the output 
of the 0AM beam splitter onto two holographic spatial light modulators (SLMs). These are liquid-crystal 
devices that can impart an arbitrary 2D phase pattern onto a photon. We use these devices to perform 
projection measurements of 0AM. These measurements are carried out by flattening the phase of incident 
photons such that they couple into a single-mode optical fiber.l^ The scheme is based on the fact that a 
single-mode fiber (SMF) only carries spatial modes with £ = 0. In order to detect a photon with ^ = 1, 
a hologram with £ = —1 displayed on the SLM projects the photon into mode £ = 0. The photon then 
couples efficiently to an SMF and is guided to an single photon avalanche photodiode (Excelitas SPCM- 
AQRH-14-FC). The same procedure also applies to projective measurements of mode superpositions. Note 
that in our experiment, phase-only holograms are used, which take only the vortex phase structure of 
a Laguerre-Gaussian mode into account, and not its amplitude. This results in an imperfect mode overlap 
at the smfISI and leads to lowered coupling efficiencies for modes with £ = zLl. This drop in efficiency can 
be somewhat compensated by adding a mode-dependent quadratic phase to these holograms, which changes 
the effective coupling mode-waist. 


(332)-Entanglement witness. In order to prove that the state is indeed a (332)-type entangled state we 
have to prove that it cannot be decomposed into states of a smaller dimensionality structure. We thus have 
to show that it lies outside the (322)-set of states, that is the convex hull of all states that can be decomposed 
into 322 and 232 states. From the measured data we can extract the fidelity to the ideal state 

|^-) = T(|0,0,0) + |1,-1,1) + |-1,1,1)) , (9) 

which we will denote as Fexp := Tr(/)exp|^)(^|)- We thus need to compare the experimental fidelity with 


14 




the best achievable fidelity of a (322)-state, i.e. 


-F’max:= Tr(cr|«')(«'|). (10) 

(7 G (oZZ j 

If Fexp > Fmax: wc Can concludc that the experimentally certified fidelity cannot be explained by any state 
in (322) and thus the underlying state is certified to have an entangled dimensionality structure of (332). 
To calculate Fmax h is useful to observe that it is convex in the set of states, i.e. the maximum will always 
also be reachable by a pure state. Furthermore, since the set (322) is the convex hull of 322 and 232, we can 
write the following 


^max = max |(T|^)p=max[ max |(T|^)p, max |(T|^)p]. (11) 

|<I>)G(322) ^ ^ |<I>)G322 ^ ^ |<I>)G232 ^ ^ ^ 

Now for a fixed rank vector xyz^ these fidelities can be bounded by noticing that 


max |(5'|$)|2 < min[ max max |(«'|^)P 

\^)^xyz rank{Tr23\^) {^\)=x rank{Tr±3\^) {^\)=y 


max |(T|^)P]. 

rank(Tri2\^) (^|)=2; 


( 12 ) 


Now inserting (12) in (11) we can compute each of the appearing terms using a previously published 
theorem.l^ It states that if a state has a Schmidt decomposition across a cut A\A given by |T) = 


^ ( 8 ) we can compute the maximal overlap with a state of bounded rank across this par- 


tition as 


x — l 


max 

rank(Tr^\^) (<l>|)=ai 


|(^-|$)|2 = ^Af, 

i=0 


(13) 


where we assumed ordered Schmidt coefficients, i.e. > A^+i. Now all we need are the coefficients for the 
Schmidt decompositions for our target state for all three partitions. For 3|12 they are and for 2|13 

and 1|23 we get Inserting these numbers we find that the maximum overlap of the target 

state with a (322) state is given by 


F — - 

X max — „ • 


(14) 


Witness Measurements. The experimental fidelity Fexp := Tr(pexp|^)(^|) determines which measure¬ 
ments are required. The projector |T)(T| projects only onto the non-zero diagonal and off-diagonal elements 
contained in the density matrix Pexp- Additionally, for the purpose of normalization, it is necessary to 
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measure all other diagonal elements in Pexp- This results in 18 diagonal and 3 unique off-diagonal elements 
that need to be measured in order to calculate Fexp. In our experiment, we can only perform projective 
measurements with SLMs. A diagonal element is given by one single projection {ijk\p\ijk) = with 

Ct '•= 0 1 Sj=-i 0 1 S/c=o 1 containing all diagonal elements for normalization. Out of the 

six off-diagonal elements, only three are unique and need to be measured: (000|p|l-ll), (000|p|-lll) and 
(-111 IpI 1-11). Note that the last off-diagonal element is only in a two-particle superposition. Hence, it can 
be measured in the standard way that two-particle two-dimensional states are usually measured. In order 
to measure the other two off-diagonal elements with projective measurements, we decompose them into 
and cTy measurements. The real and imaginary part of each element can be written as 


^[{ijk\p\lmn) 

1 

' - (d’‘ ® 

^[{ijk\p\lmn)] 


- (4* 0 4 




0 

0 4 ’“)-( 4 ' 0 4 “ 


0 4’4 ( 15 ) 


where 4’*^ = I®) (^1 + 1^) (®l 4’*’ = ^ I®) (^1 1^) (®l- The operators are also not measurable directly 

with SLMs and are therefore rewritten using the following operators: 


V+{a,b) — 1+) (+4^6) — |a) (a| + \b) {b\ + \a) (6| + |6) (a| (16) 

P_(a,6) = I-) = 14 {a\ + |6) (&| - |4 {b\ - \b) {a\ 

V+iia, b) = l+z) {+i\^a,b) = l«) («l + 1^) (^1 - * l«) (^1 + * 1^) («l 

r-i{a, b) = |-i) (-i|(„^ft) = 14 {a\ + \b) {b\ + i |4 (&| - i\b) {a\, 

where = |4 + \b), = |4 - \b), = |4 + i\b) and = |4 - i\b). These 

superposition states can be measured with SLMs in our experiment. Thus the <t operators from Eq. [T^can 
be written in the following manner: 


4’'’ = ^(A(«4)-^-(a,4) 

a;’^ = ^(V-iia,b)-V+i{a,b)). (17) 

This leads to 64 projection measurements required for measuring each of the two off-diagonal elements 
(000|p|l-ll) and (000|p|-lll). Only 16 projection measurements are required for the off-diagonal element 
(-111 IpI 1 -11) since it involves only one dimension for photon C. Summing up these measurements as well 
as the 18 diagonal elements leads to 162 total measurements. The results of these are given in extended 
data table 1. From these measurements, the overlap between the generated state Pexp and the ideal state 
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(332)-state |^) is calculated to be Fexp = 0.801 ±0.018. The error in the overlap is calculated by propagating 
the Poissonian error in the photon-counting rates by performing a Monte Carlo simulation of the experiment. 
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Extended Data Figure 1: Detailed Experimental Setup. A TiiSapphire pulsed laser source with pulse duration of 140fs 
and repetition rate of 80 MHz centered at 808nm is focused by lens LI into a /3-Barium Borate crystal (BBO) to produce 
pump pulses at 404nm via the process of second harmonic generation. In order to avoid damage, the BBO crystal is moving 
periodically with a motor (M). Lens L2 re-collimates the 404nm pump pulses, which are separated from the 808nm laser by a 
dichroic mirror (DM1). Two OAM-entangled photon pairs are produced via type-II spontaneous parametric down-conversion 
in two periodically poled Potassium Titanyl Phosphate (ppKTP) crystals (NLl and NL2) of dimensions 1 x 2 x 1mm. The 
full-width half-maximum beam waist of the 404nm pump laser at the crystals is 240p,m with an average power of 820mW. 
To prevent gray-tracking and crystal damage, the ppKTP crystals are moved periodically up and down in steps of 200)am via 
motors (M). A second dichroic mirror (DM2) separates the entangled photon pairs from the 404nm pump laser beam. Lenses 
L5 are used to collimate the down-converted 808nm photon pairs. A 4/ imaging system L4 is used to perfectly image the 
404nm pump beam from crystal NLl to crystal NL2 to ensure that the pump beam mode is exactly the same at both crystals. 
Polarizing beam splitter (PBS) are used to deterministically separate all four photons. A moving trombone system of mirrors 
(TB) is used to ensure that photons from both crystals arrive at the OAM beam splitter at the same time. A 4/ system of 
lenses (L6) compensates the extra propagation that the photon from crystal NLl undergoes, ensuring that both photons have 
a good spatial mode overlap at the entrance to the OAM beam splitter (highlighted in green). The OAM beam splitter is 
implemented in a double Sagnac interferometer configuration that allows for stable operation over several days. Note that this 
is different from the physical implementation shown in Fig. ^ (for simplicity), but has exactly the same physical outcome. A 
piezo controlled mirror (P) is scanned every hour to optimize the alignment of the interferometer and ensure that it is sorting 
OAM modes correctly. Two 4/ systems of lenses (L6) image the output of the OAM beam splitter onto spatial light modulators 
(SLMs). Two half-wave plates (HWPs) at 45° rotate the photon polarization from vertical to horizontal in paths A and D. 
Projective measurements are performed by four SLMs in combination with aspheric lenses and single mode fibers, which guide 
the photons to single photon avalanche diodes (Det A, B, C, and T). Additionally, narrowband filters (IFl and IF2) are used 
before all four detectors to ensure spectral indistinguishability between the interfering photons. 
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Extended Data Table 1: Numerical values of the measured density matrix elements. 
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